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. The procedure to find gauge invariant variables for two-parameter nonlinear perturba- 

' tions in general relativity is considered. For each order metric perturbation, we define the 

variable which is defined by the appropriate combination with lower order metric perturba- 



00 



tions. Under the gauge transformation, this variable is transformed in the manner similar to 
the gauge transformation of the linear order metric perturbation. We confirm this up to third 
^\ ' order. This implies that gauge invariant variables for higher order metric perturbations can 

04 \ be found by using a procedure similar to that for linear order metric perturbations. We also 

derive gauge invariant combinations for the perturbation of an arbitrary physical variable, 
■ other than the spacetime metric, up to third order. 

> 

o . 

<3\ , §1. Introduction 

O 

Q ■ The perturbative approach is one of the popular techniques to investigate physi- 

cs , cal systems. In particular, this approach is powerful when the construction of exactly 

soluble models is difficult. In general relativity, exact solutions to the Einstein equa- 



o 



O ■ tion are most often too idealized to properly represent the realm of natural phenom- 

I , ena, though many exact solutions are known.^ Here the constructing perturbative 

^Jjl solutions around appropriate exact solutions is a useful approach to investigate re- 

alistic situations. Cosmological perturbation ih.eov^~^ is now the most commonly 
used technique, and perturbations of black holes and stars have been widely studied 
^ ' to obtain descriptions of the gravitational radiation emitted from them.l^^^ These 

■ recent perturbative analyses have been extended to second-order perturbations, but 

in many cases, these treatments employ an expansion in a single parameter. 

In some physical applications, it is convenient to introduce two (or more) in- 
finitesimal perturbation parameters to elucidate the physical meaning of the pertur- 
bations. One typical example is the study of perturbations of rotating stars.'^~n2J 
No exact analytic stationary axisymmetric solution describing rotating stars has yet 
been obtained, at least for reasonably interesting equations of state. To treat rotating 
stars, perturbative analyses employing the "slow rotating approximation" are com- 
monly used. In this approach, the background is a non-rotating star, i.e. spherically 
symmetric star, and two small parameters, A and e, corresponding to the pulsation 
amplitude and the rotation parameter, are introduced. The pulsation amplitude is 
given by the amplitude of the metric perturbation, and the rotation parameter is 
given by e = Q/ \J GMjB? , where Q is the uniform angular velocity, and M and R 
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are the mass and the radius of the non-rotating star, respectively. In this approach, 
the first order in e describes frame dragging effects, with the star actually remaining 
spherical, and terms describe the effects of rotation on the fluidP^ Because the 
mass-shedding limit corresponds to e ~ 1, this approximation is valid for angular 
velocity f2 much smaller than the mass-shedding limit, and this approximation has 
also been used recently in the study of the instability in rotating stars. (See the 
review paper by Stergioula^^ and references therein). This example shows that 
the two-parameter perturbation theory is interesting from the viewpoint not only 
of mathematical physics but also of its applications. There are many astrophysical 
situations that should be analyzed using multi-parameter perturbation theory. 

In spite of these efforts, classical studies in the literature have not been analyzed 
taking care of the full gauge dependence and gauge invariance of the non-linear per- 
turbation theory. For example, the delicate treatment of gauge freedom is necessary 
when we evaluate boundary conditions at the surface of the matter distribution and 
the perturbative displacement of the surface (for example, see Refs. 24)). An im- 
plicit fundamental assumption in relativistic perturbation theory is that there exists 
a parametric family of spacetimes such that the perturbative formalism is built as a 
Taylor expansion in this family around a background. The perturbations are then 
defined as the derivative terms of this series, evaluated on this background.^ To 
carry out this evaluation, we must identify the points on the background spacetime 
and those on a physical spacetime that we attempt to describe as a perturbation of 
the background spacetime. This choice of the identification map is usually called the 
gauge choice.^ The important point is that this identification is not unique, i.e. 
there is a degree of freedom involved in the choice of this identification map. This 
is the gauge freedom in the perturbation theory. Clearly, this consists of more than 
the usual assignment of coordinate labels to points of a single spacetime. Further, 
the Einstein equation does not determine this gauge freedom, and we must fix this 
gauge freedom by hand or extract the gauge invariant part of the perturbations (for 
example, see Ref. 13)). This problem does not arise when this gauge freedom is com- 
pletely fixed and if a change of the gauge is not necessary to analyze or interpret the 
physical meanings of the results. Otherwise, this problem always arises. Therefore, 
it is important to clarify the gauge transformation rules of physical variables and the 
concept of gauge invariance. 

In this paper, we present the procedure to define gauge invariant variables in 
the two-parameter nonlinear perturbation theory. Recently, Bruni and coworkers 
derived the gauge transformations and introduced the concept of gauge invariance in 
the two-parameter nonlinear spacetime perturbation theory.^ They derived explicit 
gauge transformation rules up to fourth order, i.e., including any term X'^e^ with 
k + k' < 4. We follow their ideas in this paper. Although we keep in mind the above 
mentioned practical examples, we do not make any specific assumption regarding 
the background spacetime and the physical meaning of the two-parameter family. 
Instead, we assume the existence of a procedure to determine the gauge invariant 
variables at linear order. For each order metric perturbation, we define the variable 
which is defined by the appropriate combination with lower order metric perturba- 
tions. Under the gauge transformation, this variable is transformed in the manner 
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similar to the gauge transformation of linear order metric perturbation. We confirm 
this up to third order. This implies that we can always find gauge invariant variables 
for higher order metric and matter perturbations because we have a procedure to 
determine gauge invariant variables of linear perturbations as we assumed. 

Because we make no assumption concerning the background spacetime, our pro- 
cedure is applicable to various situations. Further, we note that gauge freedom 
always exists in the perturbation of theories in which we impose general covariance. 
Therefore, our procedure is applicable not only to general relativity but also to any 
theory in which general covariance is imposed. However, we cannot treat the situa- 
tion in which the change of the differential structure arises due to the perturbations. 
We also note that the procedure developed here has already been applied to clarify 
the oscillatory behavior of a gravitating Nambu-Goto strin^^ in which it is crucial 
to distinguish the gauge freedom of the perturbations and the motion of the string. 
Through such considerations, we have already confirmed that the procedure we study 
here is applicable in a specific case. 

The organization of this paper is as follows. In ^ we present the necessary 
mathematical tools, deriving Taylor expansion formulae for two-parameter families of 
diffeomorphisms. In ^ we set up an appropriate geometrical description of the gauge 
dependence in two-parameter families of spacetimes and derive gauge transformation 
rules for the perturbations. In ^ the procedure to determine the gauge invariant 
variables of nonlinear perturbations is described. The final section, ^ is devoted to 
summary and discussions. Sections |21 and |21 consists largely of a review of the work 
of Bruni et al.,'^ which is referred to as BGS2003 in the present paper. However, 
these sections include some additional explanations that are not given in BGS2003. 
In particular, we note that the representation of the Taylor expansion given in this 
paper is simpler but equivalent to that given in BGS2003. We employ the notation 
of BGS2003 and also use the abstract index notation.'^ 

§2. Taylor expansion of the two parameter diffeomorphisms 

Perturbation theories on a manifold are usually based on a Taylor expansion on 
an extended manifold of the original manifold. Taylor expansions provide an approx- 
imation of the value of a quantity at some point in terms of its value and the values 
of its derivative, at another point. Here, a Taylor expansion of tensorial quantities 
can only be defined in terms of a mapping between tensors at different points of 
the manifold under consideration. This implies that a two-parameter perturbation 
theory on a manifold requires a Taylor expansion of such a mapping given by a 
two-parameter family of diffeomorphisms on the manifold. In this section, we review 
the Taylor expansion of two-parameter diffeomorphisms developed in BGS2003, with 
some modifications to clarify the essence of their idea. 

Given a differentiable manifold Ai, we consider a family of diffeomorphisms 
characterized by two parameters on A4, X and e: 
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As emphasized by Bruni et al.;^^' the diffeomorphisms <Px,e do not form a group in 
the form ^Ai,ei ° ^A2,e2 = ^Ai+A2,ei+e2 ^^r all Xi, ej G IR (i = 1, 2). This differs from 
the usual situation for exponential mapsP^ In the generic case, we must keep in 
mind the fact that 

^Ai,.i o ^A2,e2 / <^Ai+A2,.i+.2- (2-2) 

This means that ^^A.ei in general, cannot be decomposed into the form ^x,e = ^o,e ° 
^A,0) where both c^o.e and ^A,o are one-parameter families of diffeomorphisms. Hence, 
in the derivation of the representation of the Taylor expansion of the pull-back ^ 
of ^A,ei we cannot use the representation of the Taylor expansion of the pull-back of 
the one-parameter family of diffeomorphisms?^'^ The Taylor expansion based on 
the usual exponential maps is realized as a special case of the representation derived 
here. 

The simple algebraic properties of the coefficients of the Taylor expansion of 
for an arbitrary tensor field Q leads to their representation in terms of suitable 
Lie derivatives. We start from the formal expression of the Taylor expansion of the 
pull-back 'P'^^Q, which is given by 



fc,fc'=0 L ) A=e=0 



As the properties of the coefficients of the Taylor expansion H2-3() . we stipulate that 
the operators d/dX and d/de in the bracket {*};^^^^q in Eq. (|2-3j) are not symbolic 
notation but, rather, the usual partial differential operators on IR^. The representa- 
tion of this Taylor expansion in terms of the Lie derivatives is explicitly derived in 
Appendix We note that the Leibniz rule plays a key role in the derivation of the 
representation of the Taylor expansion H2-3|) . 

In this paper, we present the expansion of the pull-back <PX^ to order X^^ 
with /c -|- A;' = 3 in terms of suitable Lie derivatives. For this purpose, we introduce 
the following set of operators where p and q are integers, on an arbitrary 

tensor field Q: 

>C(i,o)Q:={|^^vQ}^__^> (2-4) 
Ao,i)Q:={^<^A,eQ}^__^, (2-5) 
>C(2,o)Q := {^^a,.q}^_ _^ - 4,o)Q> (2-6) 
'C(i,i)<3 := yoxofX^^] ^_ _Q ~ 2 (^(i.o)^(O'i) + '^(0,1)^(1,0)) (2-7) 
>C(o,2)Q := {^^a,.q}^_^_^ - 4,1)^' (2-8) 



-^(3,0)^ := {^^a,.q} _ _ - 3£(i,o)^(2,o)Q - 4,o)Q> (2-9) 
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^{2,i)Q ■■= { a^^^'^*^}^, _g ~ 2/:(i,o)/^(i,i)Q - 'C(o,i)/:(2,o)Q 

--^(i,o)'^(o,i)'^(i,o)Qi (2-10) 
-'^(o,i)^(i,o)^(o,i)Q) (2-11) 

^(o,3)Q {£^a,.q}^_ _^ - 3/:(o,i)/:(o,2)Q - (2-12) 

As shown in Appendix the above operators >C(p g) are hnear and satisfy the Leib- 
nitz rule, and hence they are derivative operators. Because the puh-back ^ of a dif- 
feomorphism ^ commutes with contractions and the exterior derivative P^'^ the 
operators .C(p g) also commute with any contraction and exterior derivative. There- 
fore, for each of them, there is a vector field such that 

^i(,p,q)Q '■= ^{p,q)Q (2-13) 

for g = 0, 1, 2, 3. 

Using the Lie derivative (|2-13jl , we can express the Taylor expansion (|2-3|) of the 
pull-back ^ of in terms of the Lie derivatives associated with the vector fields 
g) (|2-13p up to third order in A and e: 

^l.eQ = Q + 5(1,0) Q + ^£iio,i)Q 

+ Y {"^€{2.0) +4(1,0) }<3 

+Ae I ^£.(1^) + 2 °^?(i,o) ■^5(0,1) + 2 '^5{o,i) °^C(i,o) j" ^ 

+ |- {-^5(0,2) +4(0,1)}^ 

+ "g" I ■^€(3,0) + 3-^5(1,0) ■^€(2,0) + ■^€(1,0) J ^ 

A^e r 1 

^ I °^?(2,i) + 2-^5(1, 0) "^€(1,1) + "^?{o,i) ■^€(2,0) + °^C(i,o) °^C(o,i) °^?(i,o) / ^ 

Ae^ r 1 

^ |'^?(i,2) + 2i?5(()_;ij i^g(i_;^j + £5^^ £^((, 2) + "^^co,!) °^S(i,o) °^5(o,i) J ^ 

+ {°^?(0,3) + 3"^e{0,l)°^C(0,2) + 4(0,1)} ^3 

+0^A,6). (2-14) 

Here, we note that the definitions given in Eqs. H2-4p - (|2-12p of the derivative op- 
erators g) and the expression (|2-14p of the Taylor expansion does not include 
arbitrary parameters, while that derived in BGS2003 does. In Appendix ^ it is 
shown that the parameters in the representation derived in BGS2003 can be elim- 
inated through the replacement of the generators .^q^ without loss of generality. 
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The expressions (|2^ - (|242|) and pal]) are equivalent to those in BGS2003. As 
emphasized in BGS2003, the representation of the Taylor expansion is not unique, 
and there are several different but equivalent representations. These representations 
are reduced to Eq. H2-14p through the replacement of the generators (.^p^y Hence- 
forth, we refer to the simpler representation (|2-14|) as the ^^canonical representation" 
of the Taylor expansion of the two-parameter diffeomorphisms. Further, we denote 
this expression by e(^{p q))Q when there is a need to specify the generators C^pg)! 
explicitly. 

Next, we consider the problem how to recover the one-parameter case from the 
two-parameter case when the two parameters A and e are no longer independent, 
e.g. when e = e(A). The case in which either A or e vanishes is trivial and it 
can be recovered from the above expressions by simply setting A = or e = 0. 
Another simple case in which the two parameters A and e are linearly dependent, 
i.e., e = aA (a 7^ 0), is discussed in BGS2003, there it is shown that the Taylor 
expansion of the two-parameter case is reduced to the one-parameter case through 
the replacement of the generators C^pgy Here, we consider the more generic case in 
which the infinitesimal parameter e is given by a Taylor expansion in A : 

e = e(A) = Va„ — . (2-15) 

n=0 

Because we only consider the representation ^Ae('^(pg))Q third order, we can 
restrict our attention to the expression (|2-15|) up to third order. Substituting (|2-15|) 
into the Taylor expansion of ^Q, we obtain 

^XeQ = Q + ao£(^o,i)Q + ^^C{i,o)Q 

+ ^«o{i^C,o,.,+4o.)}^ 

+OoA ^^C(iA) + 2'^f(i,o) "^Ca),!) + 2°^?{o>i)'^^{i>o)| Q 

+ 2^0-^ {'^C(i,2) + ^■^€(0,1) -^((1,1) + '^C(i,o) °^C(o,2) + ■^€(0,1) -^((1,0) -^^(0,1) } Q 
+ 2«oA^ |^C(2,i) + 2'CC(i,o)'^C(i,i) + °^C(o,i)'^C(2,o) + ^C(i,o) ■^€(0,1) '^((i.o) } Q 

+0\X,ao), (2-16) 
where the vector fields ("[J, are defined by 

C(i,o) •= ^{1,0) + ^1^(0,1)' (2-17) 
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1 

+ "'l';(0,2)' 


(2-18) 


AO 

M2,0) 


= ^(2,0) 


+ 2aiC('i^i) + «2C('o,i) + «i^(o,2)' 


(2-19) 




- ?(1,2) 


+ «l^{0,3)' 


(2-20) 


AO 

M2,l) 


- ?(2,1) 


+ 2aiC(\,2) +«2C('o,2) +01^(0,3) + Ol[K(0,l),'^{l,0)]7C(0,l)]"7 


(2-21) 


AO 

M3,0) 


— 

- M3,0) 


+ «3C{'o,i) + 3a2C('i,i) + ^'^1^(2,1) + 3aiC{'i,2) 





+3a2ai^('Q^2) + «i'?(o,3) + 2"2K(o,i)>^(i,o)]'' 

+2ai[['^(o,i)i'^(i,o)])'^(o,i)]" + ai[K(i,o)>C(o,i)]'^(i,o)]''- (2-22) 



Equation (|2-1H|) has the form of a Taylor expansion of a diffeomorphism with two 
infinitesimal parameters, A and gq. Equation H2-16() shows that the coefficient ao in 
the expansion ()2-15|) plays the role of an infinitesimal perturbation parameter that 
is independent of A. When ao = 0, Eq. H2-16|) reduces to the Taylor expansion in 
the case of a single infinitesimal parameter A. Thus, even when two infinitesimal 
parameters depend on each other as in the relation H2-15() . we find that the Taylor 
expansion H2-14() is reduced to that of a single parameter, as in the trivial case e = 
or A = 0. 

We next derive the representation of the inverse of the canonical representation 
<j))^ arbitrary tensor field Q. To do this, we first consider the product 

'^A.e ° ^A,e of the two diffcomorphisms \^x^i: and ^^A.e- We consider the canonical 
representations ^XeiC(^pq))Q ^'^d '^*x e(^{p q))Q- obtain the Taylor expansion of 
the pull-back of °^\,e: we first derive (p*^ ei^fp g))*^ arbitrary tensor field S 

and substitute the canonical representation of the Taylor expansion 5 = 'I^x^iC^p q))Q- 
Then, we obtain the representation of the pull-back, 

{^xAC^p,q)) ° '^xAiM)yQ = ^U^m) ° n.iCM^Q- (2-23) 

In order that we can regard '^^XeiC^pq)) inverse of ^Xei^^pq))' '^^ choose the 

generators C(p^g) so that (^^\,e{C{p,q)) ° '^x,e{^(p,q))) Q = Q- This is accomplished by 
choosing 



Cfi,o) - 


-?(1,0)' 




(2-24) 


Cfo,i) = 


~MO,i)' 




(2-25) 


km - 


'^(2,0)' 




(2-26) 


Q"o,2) = 


ca 
~?(0,2)' 




(2-27) 


Afl _ 


pa 




(2-28) 


= 


-^(3,0) + 3[C(2,0),C(1,0)]", 




(2-29) 


C{2,1) = 


-'^"2,1) + K(2,0)i'^(0,l)]" + 2['^(1,1) 


^(1,0)]") 


(2-30) 


fa _ 


"■^^,2) + K(0,2)i'^(l,0)]" + 2^(1,1) 


^(0,1)]''' 


(2-31) 


Q'o.s) = 


-^(0,3) + 3[C(0,2),C(0,1)]''- 




(2-32) 
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Then, the exphcit form of (<^A,e(?('p.g)))*<3A,e is given by 
= Q- Ai:^(^ - e£(.^^^^Q 

+y {--^5(2,0) + 4(1,0)} '3a,. + Y {--^€(0,2) + 4(0,1)} '3 

+Ae | — + 2°^?(i,o) °^?(o,i) ~^ '^^i(OA)^i{i,o)^Q 



6 



?(3,0) '^■^5(2,0) °^5(1,0) '^^( 



a,o)}^ 



Ae2 



£<. 



2 I -€(2,1) °^ 5(2,0) "^€(0,1) ^°^C(i,i) ■^€(1,0) °^€(i,o)°^?{ 



+ - 



{ 



■^€(1,2) + °^5(a,2) °^?(i,o) + 



£pii~. -,^ £pi^ 



5(1,1) €(0,1) °^5(0,1)°^?(1,0) C(i 



0,1)°^5(1,0) } ^3 

.0,1) }^? 



+ {"■^€(0,3) + 3=^5(0,2) -^€(0,1) - 4(0,1) } ^ + ^) 



(2-33) 



This exphcitly shows that ^^^j; 7^ ^_A,-e, as emphasized in BGS2003. Further, if all 

generators i^^^^^ commute, we obtain the equality (^A,e)*(^"p,g))Q = ^XS.~%,q))^- 
Finally, we show that the two-parameter group of diffeomorphisms that satisfy 
the property 

0Ai,£i o 4>\2,e2 = 0Ai+A2,ei+e2 VA, e G IR (2-34) 

is obtained as the special case of the two-parameter family of diffeomorphisms. This 
property implies that the two-parameter group (/)A,e can be decomposed into two 
one-parameter groups 4>\ o aiid e of diffeomorphisms: 



9A,e = (PA.O ' 



90,£ o <Pa,o- 



(2-35) 



These two one-parameter groups of diffeomorphisms are generated by the vector 
fields r/p^^ and r/^^^, respectively. Each of these vector fields is defined by the action 
of the corresponding pull-back, 0^ q and (p^ for a generic tensor field Q on x M: 

:= lim y(<aI,oQ - Q). £v,.,Q ■= li™ - Q)- (2-36) 

Because the property (|2-34|) implies that the two-parameter group is Abelian, 
the vector field 77"^^ and ry"^^ commute 

M,%)]" = 0. (2-37) 

The Taylor expansions of the pull-backs (f)\ qT, 0q are given 



00 ^fc 



fc=0 



fe=0 



00 



A=0 fc=o 



-^=0 fe=0 



(2-38) 
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Then, using the decomposition (|2-35|) . we obtain the Taylor expansion of the two- 
parameter group of pull-backs 0^ ^Q: 



k^k' 



k,k'=0 



klk'l 



Qk+k' 
QXkQfk 



7<PleQ 



A=£=0 k,k'=0 



This expression is also obtained as the special case of the Taylor expansion (|2-14() 
of the two-parameter family of diffeomorphisms ^\ imposing the conditions that 
^(A) = ^(1,0) ^ 0' ^(e) = ^(0,1) / 0, [^(1,0), e(o,i)]" = Oi and i^^^^ = for p + g > 1. 



§3. Gauge transformation of perturbation variables 



Using the above Taylor expansion of two-parameter diffeomorphisms, we con- 
sider gauge transformations in two-parameter perturbations of the manifold. 

3.1. Gauges in perturbation theory 

Let us consider the spacetime (A^O) ^^^9ab), which is the background spacetime for 
the perturbations, and a physical spacetime {M.,gab), which we attempt to describe 
as a perturbation of the background spacetime {A4o, '^^^Qah)- Let us formally denote 
the spacetime metric and the other physical tensor fields on the physical spacetime 
M by Q. In perturbation theory, we are used to write expressions of the form 

Q{x) = Qq{x)+5Q{x). (3-1) 

This expression relates the variable Q on to the background value of the same 
field, Qoj and the perturbation 5Q. In the expression H3-l() . we have implicitly 
assigned a correspondence between points of the perturbed and the background 
spacetime. This is the implicit assumption of the existence of a map Mq M : 
p € A^.^ This correspondence associated with the map TWq -M. 

is what is usually called a "gauge choice" in the context of perturbation theory.*) 
Clearly, this is more than the usual assignment of coordinate labels to points on 
the single spacetime. It is important to note that the correspondence established 
by such a relation as Eq. (|3-1|) is not unique. Rather, Eq. (|3-1|) involves the degree 
of freedom corresponding to the choice of the map Mq — > M. (the choice of the 
point identification map A^o ~^ M)- This is called "gauge freedom". Further, such 
freedom always exists in the perturbation of a theory in which we impose general 
covariance. 

Here, we introduce an (m + 2)-dimensional manifold M to study two-parameter 
perturbation theory based on the above idea. The manifold M is foliated into m- 
dimensional submanifolds diffeomorphic to Al, so that Af = A4 x IR^. Each copy of 
A4 is labeled by the corresponding value of the parameters (A, e) € M^. The manifold 

*•* More precisely, as mentioned in BGS2003, Eq. 13-11 1 gives a relation between the images of the 
fields in IR™ rather than between the fields themselves on the respective manifolds M and A4o, i.e., 
we are saying that there is a unique point s £ IR that is at the same time the image of two points: 
one in Mo and another in Ml. However, in this paper, we deal with only the point identification 
map X : Mo M. 
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J\f has a natural differentiable structure that is the direct product of that of A4 and 
]R^. With this construction, the perturbed spacetimes Mx^e for each (A, e) must have 
the same differential structure, and the changes of the differential structure resulting 
from the perturbation, for example the formation of singularities, is excluded from 
our consideration. Each point on M is assigned by (p, A,e), where p E A^A,ei and 
each point on the background spacetime Mo in M is assigned by A = e = 0. 
Let us consider the set of field equations 

£[Qx,e] = (3-2) 

on TWa.e for the physical variables Q on A^A,e- The field equation H3-2() formally 
represents the Einstein equation for the metric on Ai\^e and the field equations 
for matter fields on A^A,e- If a tensor field Qa.c is given on each A^A,ei Qx,e is 
automatically extended to a tensor field on TV by Q{p, A, e) := Qx^eip)-, with p G M.\^^. 
In this extension, the field equation (|3-2j) is regarded as the equation on TV. 

Now, we define the perturbation for an arbitrary tensor field Q by comparing 
Qa.e with Qq. To do this, it is necessary to identify the points of ^A\^^ with those of 
Mq. This is easily accomplished by assigning a diffeomorphism AfA,e : M ^ M such 
that X\^e ■ Mq —>■ TWA,e- It is natural to regard Xx^e as one of the two-parameter 
group of diffeomorphisms that satisfy the property ()2-34|l . Then, is generated 
by two vector fields "^V^x) and '^^/('j,) on TV that satisfy 

[''^(A),%)r = 0. (3-3) 

Further, the normal forms of A^A.e in TV are given by {dX)a and {de)a, and their 
duals are defined by 

|;)°(<iA), = i, = (|;)°(*). = «. (D^"^)"""- 

(3-4) 



The vector fields "^^^x) ^^'^ ^^(e) chosen so that 
where and 0"^^ are tangent to TV4A,e for each A and e : 

^{A),(.)(^e)a = ^(A),(.)(^A), = 0. (3-6) 

The choice of the vector fields 9^^-^ is essentially arbitrary, except for the conditions 
1)3 -31) and H3-6p . Therefore, we choose 0"^^ so that 

^^^tx)M = £§-/hw = (3-7) 

for simplicity. Except for these conditions, we can regard ^"^^ ^^^^ as essentially arbi- 
trary vector fields on TWa^e (not on TV) that all commute. 
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The perturbation AqQx^^ of a tensor field Q for a gauge choice <Y can now be 
defined simply as 

A^Qx^,:=A:1,Q\^^-Qo. (3-8) 
The first term on the right-hand side of (|3-8|) can be Taylor-expanded as 



oo 



klk'l ^ 

k,k'=0 



where 



^ieQU= E ^^r^Q, (3-9) 

Qk+k' 



Mo 



(3-10) 



J A=e=0 

Equations ()3-8() - H3-10p define the perturbation of order (fe, k') of a physical variable 
Q for the gauge choice X and its background value 8^^''^^Q = Qq. 

3.2. Gauge invariance and gauge transformations 

Let us now suppose that two gauges X and y are generated by the pairs of vector 
fields ('^?7(A)) '^?7{e)) and (''^??{A)i'^^(e))) respectively. These vector field are defined on 
M as the vector fields in Eqs. H3-5|) : 

/ f) \ "' / f) \ "- 



The condition ()3-3() for each set of generators, "^^^p^) and ~^'^'^^X) (e)' i^pli^s that 
the two-dimensional tangent space spanned by ^'^'(X) (e) ('^^(a) (e)) possesses a two- 
dimensional integral surface. These integral surfaces of (^-j and define 
two two-parameter groups of diffeomorphisms X and y on H . Further, ^^-j and 

'^^(A) (e) everywhere transverse to A^A,ei and points lying on the same integral 
surface of either of the two are to be regarded as the same 'point within the respective 
gauge. (See Fig. 0) Then, X and y are both point identification maps. When 
^(A) (e) 7^ ''^A) (e)' t^ese point identification maps are regarded as two different gauge 
choices. 

The pairs of vector fields "^rf^ys) (e) "^^(A) (e) ^'■^ both generators of two- 
parameter groups of diffeomorphisms, and they pull back a generic tensor Q to 
two other tensor fields, X'^J^ and y\ ^Q, for any given value of (A, e). In particular, 
on A^Oi we now have three tensor fields, i.e. Qq and the following two : 

""Qx,. ■■= KMmo ' ^^^.^ ■= ^IMmo ■ (^-i^) 

Because X and y represent gauge choices mapping the background spacetime TWq 
into a perturbed manifold A^A.ei as mentioned above, "^Qx^e and ^Q\^e are the rep- 
resentations in Alo of the perturbed tensor for the two gauges. Using (|3-8|) - (|3-1U() . 
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Fig. 1. Point identification from tlie background manifoid to tlie perturbed manifold. 



we can write 



k,k'=0 



k,k'=0 



(3-14) 
(3-15) 



where S^p^^Q, and dy'^'^Q are the perturbations (|3-in|) in the gauges X and 3^, 
respectively ; 



Jk,k') 



Mo 



Ak,k')^ _ nk nk' ^ 



Mo 



(3-16) 



Following Bruni et al. we consider the concept of gauge invariance up to order 
(n,n'). We say that Q is gauge invariant up to order (n,n') iff for any two gauges 
X and y 

d'^-^'^'^ Q = Sf^'^ Q y{k,k'), with k<n, k'<n'. (3-17) 

From this definition, we can prove that the (n, n')-order perturbation of a tensor field 
Q is gauge invariant up to order (n, n') iff in a given gauge X we have £^'^5^'^'^ = 
for any vector field defined on A^o and for any (k, k') < (n, n'). As a consequence, 
the (n, n')-order perturbation of a tensor field Q is gauge invariant up to order (n, n') 
iff Qo and all its perturbations of lower than (n, n') order are, in any gauge, either 
vanishing or constant scalars, or a combination of Kronecker deltas with constant 
coefficients.E2»JIlM» 
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Next, we consider the gauge transformation of a tensor field Q. If a tensor Q 
is not gauge invariant, its representation on Mq does change under a gauge trans- 
formation. To consider the change under a gauge transformation, we introduce the 
diffeomorphism <Px,e '■ -Mq — > TWq for each value of (A, e) E IR^. The diffeomorphism 
^A,e is defined by 

^A,. := {Xx,ey' O J^A.e = ^-X,~e o 3^A,., (3-18) 

where we have used the fact that the point identification map X is a two-parameter 
group of diffeomorphism. When X and 3^ are regarded as two different gauge choices, 
<^A,e represents the gauge transformation from the gauge X to the gauge 3^. It is 
important to note that as a consequence of the definition H3-18() . ^A,e '■ -^o ~>' 
Aio does not become a two-parameter group of diffeomorphisms in A^Oi while the 
identification maps X and y are both two-parameter group of diffeomorphisms. 
Actually, ?^Ai,ei °^\2,e2 7^ ^Ai+A2,ei+e25 "^^^ t° ^hc fact that the vector fields '^r?p^) j-^) 
and '^f]'(x) (g) have, in general, a non-vanishing commutator. In spite that ^A,e is not a 
two-parameter group of diffeomorphisms, the representation of the Taylor expansion 
of the pull-back for an arbitrary tensor field Q is obtained by using the results 

of 

The tensor fields '^Qx,e and '^Qx.e, which are defined on Mo by the gauges X 
and y, are connected by the linear map <P\ ^ as 

= ^v('^'A,eQ)Uo='^A/QA,.. (3-19) 

Therefore, the gauge transformation to an arbitrary order (n, n') is given by the 
Taylor expansion of the pull-back ^Q, whose terms are explicitly given in ^ Up 
to third order, the explicit form of the Taylor expansion is given by 

+y {"^€{2,0) + 4(1,0)}^ 

+ |-{-^C(0,2) +4(0,1)}^ 

^3 , ^ 

+ Y |°^?{3,0) + 3"^«(1,0) °^?{2,0) + 4(1,0)/ '5 

A^e r "1 

"'""^ |'^?(2,l) "I" ^°^C(1,0) + ■^C{0,1) °^?(2,0) "I" °^?(1,0) '^?(0,1) °^?(1,0) J ^ 

Ae^ r 1 

^ ^ I °^?(l,2) + ^°^C(0,1) ~'~ °^5{1,0) ■^€(0.2) + "^^CO.l) °^?(1,0) °^5(0,1) / ^ 

+ q" {4o,3) + 3"^«(o,1)"^?(0,2) + 4(0,1)} ^5 

+0\X,e), (3-20) 



14 



K. Nakamura 



using the canonical representation (|2-14|) of ^Ae(^(p(j))' where C(p,g) are now the 
generators of the gauge transformation ^A,e- 

Comparing the representation (|3-2flj) of the expansion in terms of the generators 
^{p,g) °^ pull-back and that in terms of the generators "^ilfx) [e) ^^'^ '^^(A) (e) 

of the pull-back 3^^ e ° '^-x ~eQ (~ e^)' ^^'^ explicit expressions for the gen- 
erators 1^1^ of the gauge transformation $ = o y in terms of the gauge vector 
fields '^rj^x) (e) ^^'^ '^^(x) (e) ■ Here, we give their expressions up to third order: 



ca _ y a _ X a 
?(1,0) - ^(A) ^(A) 
,a na 

- \X)-^(XY 

to. _ y„a X a 
?(0,1) - \e) - %) 
,a na 

e(2,0)=[^{A),^{A)]'^ 



[^(A), 



e(l,l) = ^[^(A)/%)]" + ^[%),^^{A)]'^ 

= 2^^{X),i^(e)T + ■^[G{X),l^(e)T-, 
e(0,2)=[%)/^(e)]'^ 

er3,0) = [[''^{A),^(A)],^(A)-2^77(A)r 

= [[^(A), '-(A)], '-(A) - 26'(A)]'', 

^(2,1) = [[^(A)/%)] /.(A)]^ [[^(A),^(.)] ,^(A) -^(A)]" 

= [ V(X) , ^(e)] , ^(A)] " + [ [^(A) , ^(c)] , '•(A) - ^(A)] " , 

= [[^(.),%)] ,^.(A)]V ,^(.) -^(.)]'' 

= [ [He) > ^(e)] > ^(A)] " + [ [0{e) , ''{.)] , ^(c) " ^(e)] " , 

e(0,3) = [[%)>^{.)]/%)-2%)]'' 



3-21) 
3-22) 
3-23) 

3-24) 

3-25) 
3-26) 
3-27) 
3-28) 
3-29) 



The expression ()3-24|) of the generator ^"-^ is different from that derived in BGS2003. 
This is due to the difference between the representations of the Taylor expansion of 
the pull-back <P*^ ^Q. In the perturbation theory, these expressions, H3-2ip - (|3-29p . are 
evaluated on the background spacetime A^o- Then, these expressions show explicitly 
that the generators of the gauge transformation (Px,e = <Y~^ o3^ are vector fields 
on the background A4q. Further, Eqs. H3-23p - ()3-29p show that the generators Cq,^) 
with p + q > 1 naturally arise from the non-commutativity of the gauge generators 

^Vlxue) and ^Vlx),ie)- 

We can now derive the relation between the perturbations in the two different 
gauges. Up to order (n, n') with n + n' < 3, these relations can be derived by 
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substituting (13^11) and (|3T5|l into 



4''°^Q-4''°^Q = %,o)Qo, (3-30) 
4°''^Q-4°''^Q = i?C(o,i)Qo, (3-31) 
4^'°)q - 4^'°)q = 2£5^^ .,,5j'°)Q + + 4, „J Qo, (3-32) 



+ |°^5(i,i) 2°^?(i.o)'^5(o,i) + 2 "^^(o^i) "^^(I'O) } (3-33) 
s'fi'^)Q - 40.2)q ^ 2£5^„^,4°'^)q + + go, (3-34) 



+ {°^C(3,o) + 3"^€(i,o) ■^e{2,o) + -^1(1,0) } Qo, (3-35) 



4^'^)q - = 2=^5(1,0) + -f?(o,i)4''°^g 



+2|"^C(i,i) + 2°^?(i,o)°^?(o.i) + 2°^€{o,i)°^5(i,o) I "^i"' ''^ 

+ ■^€(1,0) -^^CO,!) °^?(l,0) }Qo, (3-36) 

+2 + 2°^^(i'"''^^(°'i) ~'~ 2°^^("'i)°^^(i''') } ""^ 

+ °^C(o,i) "^5(1,0) "^€(0,1) } '3o, (3-37) 

+ { £^(0.3) + 3-^?(o,i) °^C(o,2) + °^|(o,i) } '3o- (3-38) 

These results are, of course, consistent with the gauge invariance up to order (n, n') as 
introduced above. Equation (|3-3Ujl [or (|3-31|) ] impHes that Qa.c is gauge invariant up 
to order (1, 0) [or (0, 1)] iff £^Qo = for any vector field on Mq. Equation (HT^ 

implies that Qa.e is gauge invariant up to order (2, 0) iff £^Qo = and £^6^^'^^Qo = 
for any vector field on Mq. This can be repreated analogously for all the orders. 
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3.3. Coordinate transformations 

The above formulation of the perturbations and their gauge transformation are 
independent of the expHcit form of the coordinate system. In some situations, it is 
convenient to introduce an exphcit coordinate system in order to carry out exphcit 
calculations in a practical case. Then, it is instructive to show the above geometri- 
cal formulation of gauge transformations in terms of the corresponding coordinate 
transformations, though it is not necessary for the discussion in §§4 and 5. Here, we 
briefly discuss this coordinate transformation. The explicit forms of this transfor- 
mation given below show the fact that the gauge freedom in perturbation theory is 
more than the usual assignment of coordinate labels to points on the single space- 
time. Details of the outline given here are given in the series of papers by Bruni and 
coworkers ^ 

We have considered two gauge choices, represented by the groups of diffeomor- 
phisms and yx^e, under which a point o on the physical spacetime A^A,e corre- 
sponds to two different points in the background manifold Mo- P = -^xli^) and q = 
y^lio) as depicted in Fig. ^ The transformation relating these two gauge choices 

is described by the two-parameter family of diffeomorphisms ^A,e = ° y\,e, so 
that = P- Under this gauge transformation, a tensor field Q on p € Mq is 

pulled back to the tensor field {<P*Q){q) = fP*{Q{p)) on q e Mo- 

Now, let us consider a chart {U, X) on an open subset U of the background Mq. 
The coordinate system X is a map from the manifold A^o to IR*". Since the gauges 
Xx^^ and 3^A.e are maps from the background Mo to the physical spacetime M\^ei 
these gauges define two maps from Mx^e to IR™: 

X o : Mx,e ^ K'" Xo y-^^ : Mx,e ^ K'" 

oi— >x(p(o)), o I—I- x{q{o)). (3-39) 

The gauge transformation 'Px,t is regarded as the transformation of these maps. 

It is well-known that there are two different points of view with which we can 
regard the gauge transformation (Px.t as a change of the coordinate system, the active 
point of view and the passive point of view. In the active point of view, one considers 
a diffeomorphism that changes the point on the background Mo, and the coordinate 
change x^{p) x^{p) is given according to the definition of the pull-back of x as 

x^'{p) := xf'i^ip)) (3-40) 

By contrast, in the passive point of view, we introduce a new chart {W, Y) defined 
hy Y := X o and the two sets of coordinates are related by 

y^ip) = x^{q). (3-41) 

In this passive point of view, the gauge transformation is regarded as not changing the 
point on A^o but changing the chart from (U, X) to {W , Y) (i.e. changing the labels 
of the points on A^o)- The coordinate transformation is given by x'^{q) y^{q)- 

Now, let us consider the transformation of a vector field V"" and the coordinate 
transformation from the active and passive points of view. 
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From the active point of view, the components of a vector field V"" in the 
chart {U, X) are related to the components of the transformed vector field V"" as 

= (X.V)* = (X.^l^V)". (3-42) 

In order to write down explicit expressions, we apply the expansion of the pull-back 
of ^ [see Eq. (|2-14() ] to the coordinate functions x^. Then, the active coordinate 
transformation is given by 

y \H0,2) "^''(0,1)^''^(0,1) 



0) 



Ae^ 



£3 



+0'(A,e), (3-43) 

where the vector fields and their derivatives are evaluated in x{p). This ex- 

pression gives the relation between the coordinates, in the chart (p{,X), of the two 
different points p and q oi Mq. 

From the passive point of view, we can use the properties relating the pull-back 
and push- forward maps associated with diffeomorphisms, 

X,^l,V = X,<P;l^V = Y,V. (3-44) 

Thus, we obtain the well-known result that the components of the transformed vector 
in the coordinate system X are defined in terms of the components of the vector 
V"" in the new coordinate system Y: 



V^^{x{p)) = {Y.V{q)r = V'^^iviq)) ' 



dx'' 



V'ixiq)). (3-45) 

x{q) 



In order to write down explicit expressions, we apply the expansion of the pull- 
back of (p'^^ [see Eq. (j2-14|) ] to the coordinate functions x^. The passive coordinate 
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A 
+— 



2 



transformation is found by inverting (|3-43|) . Using the representation (|2-33|) of the 
Tayfor expansion of the inverse of we obtain the passive coordinate transfor- 
mation: 

yf^iq) :=x^(p) x^) (q) 

= x^(<7)-Aefi_o)(<7)-eefo,i)(^) 

{-er2,o)(5)+^(i,o)(^)9.ef;,o)(?)} 

{-ero,2)(5)+^(o,i)('?)^'^^ro,i)(9)} 

+y{-^(3,0)(^)+ 3^(2.0) (?)9.e[;,o)(5) 

-efi,o)('?)5p (^(1,0) (9)^<i,o) ('?))} 

+^ {-^f2,i)('?) + ^(2,o)(9)5.e[;,i)(9) + 2efi,i)(g)9.ef^ 

-^(i,o)(^)^p(e(o.i)('?)^-^fi,o)(^))} 

+ ^ {-^il,2)i<l) + Cro,2)(9)5.efi,o)('?) + 2^(1,1) ('/)5^efo,l) (9) 

-^(o.i)(9)^p(e(i,o)('?)s.^ro,i)(^))} 

-efo,i)('?)5p (^(0,1) (9)5.^(0,1) ('z))} 

+0^(A,e). (3-46) 

This gives the relation between the coordinates of any arbitrary point q E Aio in the 
two different charts {U,X) and {U',Y). 

§4. Gauge invariant variables of higher order perturbations 

Now, we consider the definitions of gauge invariant variables for the perturba- 
tions. The gauge invariance we consider here is that up to order (n, n') for (n, n')- 
order perturbations as mentioned in ^3.21 We do this because the gauge invariance 
to all orders is not so useful. Of course, the definition of gauge invariant 
variables is not unique, because any function of gauge invariant variables is also 
gauge invariant. In this section, we present one procedure to define gauge invariant 
variables. To do this, we first consider the procedure to obtain the gauge invariant 
variables for metric perturbations. Next, we extend the procedure to define gauge 
invariant variables for any physical variables, other than the metric. 
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4.1. Metric perturbations 

As seen in Eqs. (|3-14j) and (|3-15j) . we first consider the Taylor expansion of the 
spacetime metric ga6(A, e) on Aix^e- As discuss in ^3.21 the Taylor expansion of the 
spacetime metric up to third order is carried out by using a gauge choice X on the 
background spacetime A^o^ 



k,k'=l 



gab - 9ab^ x'^ab 



— 9ab + A-^-afe + *^ + ^''''af' + x'^a-b + 

(3,0) A^e (2,1) Ae2 (1,2) (0,3) 

+0(A,e)^ (4-1) 

From Eqs. H3-3Up - ()3-38() . the gauge transformation rules for the metric pertur- 
bations up to third order are given by 

^''>afe-^''5/la6 = %,o)^°^ffafe, (4-2) 
(0,1), (0,1) , _ n (0)„ ox 

^''%b - ^'^hab = 2£^^J'^'Ik, + {^5^, .,, + £l, ,]^'^9ab, (4-4) 

(1,1), _(1'1)/, _ (0,1) L I p (1,0), 

+ + 2°^€{i.o)°^«(o,i) + 2°^«(o,i)'^«(i.o)} ^^^S'^'" (^■^) 

= 2£^^^J'^'Ik, + {i:^,^^, + £\,,]^'^9ab, (4-6) 
(3,0) _(3,0), (2,0), ,o\p , p2 \(1,0)l 



(0,2), _(0,2), (0,1) 
y'^ab X <^b 



yiT'ab 



(2,1), _(2,1), (1,1), ,n (2,0), 



+ {%,o) + 4(1,0) 



1 „ „ 1 



i,o)}^°^5a., (4-8) 



(1,2), _(1,2), (1,1), n (0,2), 



+ {^«(0,2) +4(0,1)} 



+^ '^5{l,i) + 2'^^(i'0)'^^(0^i) 2°^^('''i)°^^(i'0) \ ~ x"'ab 



(0,1) , 
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"'"'^?(0,l)'^?(i,0)'^€{0,i)} ^'^''dab, (4-9) 



(0,3), (0,3), (0,2), ,o(p , p2 \{0,1)> 

+ {^€(0,3) + 3i?5(o,.) ^5(0,2) + ^lo,)} ^^-^^^ 

Inspecting these transformation rules, we consider the procedure to separate the 
gauge invariant parts and gauge variant parts of the metric perturbation at each 
order. The aim of this paper is to show that this separation for higher order per- 
turbations can be carried out with the same procedure as for hnear perturbation 
theory. In other words, if we are able to carry out this separation at linear order, 
the separation for higher order perturbations can also be carried out. 

4.1.1. Linear order perturbations 

Suppose that, inspecting the gauge transformation rules H4-2|) and ()4-3() . the 
0(A) and 0(e) perturbations can be decomposed as 

(1,0), _. (1,0)^ (1,0) (1,0)^ 

(0,1), _. (0,1)^ (0,1) (0,1)^ 

SO that the variables ^^'^7^a6 and ^^'^^Hab are gauge invariant; i.e., 

(1,0)^ (1,0)^ _^ (0,1)^ (0,1)^ _^ 

under the gauge transformation <P = X^^ o y. Next, note that the vector fields 
and a are transformed as 

(1,0) (1,0) _ ca (0,1) (0,1) - ■ (A 

yX - -^A -^(i,o)> - x^ ~MO,i)' '^^''-■^) 

i.e., ^^'^X'^ and ^^'^^X" are the gauge variant parts of the metric perturbations 

^^'^xhab and '''^'^/lab, respectively. We also note that the number of independent 

components of ^"^''^T^ab (or ^^'^x'^a.b) is m{m— l)/2, where m is the dimension of the 
spacetime manifold. 

It is non-trivial to carry out the systematic decomposition ()4-ll|) on an arbitrary 
background spacetime, and the procedure completely depends on the background 
spacetime (Mo^^^^ dab)- simple background spacetimes in which there are some 
Killing symmetries, one useful type of analyses to carry out the decomposition ()4-ll|) 
is that based on the expansion in harmonic functions on a submanifold of the back- 
ground spacetime {Mo,^^^ gab)^^^^ For example, the harmonic functions on a 
homogeneous and isotropic three-dimensional space are used in cosmological pertur- 
bation theory.EJ Analyses based on a harmonic expansion depend strongly on not 
only local symmetries of the background spacetime but also the global topology of 
the submanifold on which scalar, vector, and tensor harmonics are defined. 
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In spite of the non-trivial nature of this decomposition, we assume that the de- 
composition H4-lip can always be carried out in some manner. What is necessary for 
our discussion here is only the result of the extraction of m gauge variant components 
from m(m + l)/2 components of the metric perturbations ^^'^/iafe or ^'^'^/lafe- The de- 
tails to realize the decomposition (|4-11() are not important to our discussion, though 
the existence of the procedure is crucial. As seen below, if such a procedure exists, 
we can easily show that the procedure to define the gauge invariant variables of 
higher order metric perturbations is reduced to that for linear perturbations, whose 
existence we assume. Though a generic formula to define gauge invariant variables 
for any order perturbation might exist, we give only the formulae for two-parameter 
perturbations up to third order. 

4.1.2. Second order perturbations 

Here, we consider the definitions of gauge invariant variables for O(A^), O(e^) 
and 0(Ae) metric perturbations. 

First, we consider the O(A^) metric perturbation. The metric perturbation 
^'^'^^hab of this order is transformed as in Eq. 1)4 -41) under the gauge transforma- 
tion <^ = o 3^. Inspecting this transformation rule, we define a variable ^'^'^^TCab 
with the gauge X by 

^''°i^afe := ^''°i/ia6 - 2^(1,0)^^ + -C^o^^ ^°^5a6. (4-14) 

We also define a variable ^'^'^yHab with the gauge 3^ by simply replacing X with y 
in Eq. (|4-14() . Using the gauge transformation rules H4-2() . (|4-4p . and ()4-13() . we can 

easily check that the variable ^'^''y'Hab is transformed as 

y'Hab ^ x'^ah ~ ^<y(2fi)^ ^ 9ahi (^'15) 

where the vector field tr^g o) defined by 

^(2,o):=er2,o) + K(i,o)/''°]^]'^. (4-16) 

We note that the gauge transformation rule (|4-15p has the same form as the gauge 
transformation rules ()4-2() and (|4-3() of the linear metric perturbations. Because 
we assume that the procedure to decompose a tensor field of second rank, which 
transforms as ()4-2() or ^A-'6\ . into the form ()4-ll() exists, the transformation rule 

(|4-15p implies that we can decompose the tensor ^^'^T^ab into a tensor field ^"^'^^Hab 
and a vector field '•^'^X* as 

^'''in^b =: ^'''in., + vJ^'^iXb + W'^ix^, (4-17) 
where the tensor ^"^'^^Hab is gauge invariant, i.e., 

^'^'^y'Hab - ^'^'^x'^ab = 0, (4-18) 
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and the vector ^^''^Xa is gauge variant, i.e., 

- (^'^X'^ = ^^2,0) + [^{1,0), ^'ixr. (4-19) 

Thus, we can extract the gauge invariant part ^"^'^^Hab of the O(A^) metric pertur- 
bation using the procedure to define the gauge invariant variables for hnear order 
metric perturbations. 

The extraction of the gauge invariant part from the O(e^) metric perturbation 
is accomphshed in a manner that is completely parallel to the above procedure for 
O(A^) perturbations. First, we define the variable '•'''^T^afe by 

^°'5^afe := ^°'5/la6 - 2^(0,1) ^^°''>a6 + £^.1) ^^^^Qab, (4-20) 

with the gauge X. From the gauge transformation rules (|4-3() . H4-6() and H4-13p . the 
variable ^^''^^Tiab is transformed as 

(0,2)^ _{0,2)^ (0)^ (.r,-,^ 

y I'-ab X "-^ — ""(0.2) ilahi ^-LJ 

under the gauge transformation A' — > 3^, where 

^(o,2):=ero,2) + K(o,i)'^°'']^]''- (4-22) 

Because the gauge transformation rule (|4-21|) has the same form as that for the 
linear-order metric perturbations, H4-2() and H4-3() . we can decompose the variable 

^ ' x'Hab as 

(0,2)^ _. (0,2)^ ^ (0,2) Y (0'2) y (A 

where ^^''^^Hab is gauge invariant, and the vector field '•'^'^Xq is the gauge variant 
part of the 0{e'^) metric perturbations. The vector field ^^''^■^Xa is transformed as 

("'JX'^ - ("'^X'^ = Cf,^,^ + [e(o,i)/°''ix]'^ (4-24) 

under the gauge transformation X ^ y. Thus, we can find the gauge invariant part 
^^'"^^"Hab of the O(e^) metric perturbation. 

In a similar way, we can also define the gauge invariant variables of the 0{eX) 

metric perturbation ^^'^/lat- The O(eA) metric perturbation is transformed as in 
Eq. H4-5|) under the gauge transformation X ^ y. Inspecting this gauge transfor- 
mation rule, we first define the variable ^^'^^Ti.ab by 

(1,1) -p .-(1.1);, f (1-0), n (0,1). 

x^ 

+ 9 {"^(1.0) ^(^"^(0,1) + £{0,i) ^■£(1,0) ^1 ^^^gab (4-25) 

in the gauge X . The variable ^^'y'Hafe is defined in the gauge y in the same form as 
Eq. H4-25() . with the replacement X ^y . Using the gauge transformation rules ()4-2p . 
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(|4-l-{)) . (|4-5|) and (|4-13j) . we can easily check that the variable ^^'^^Tlab is transformed 
as 

^^'^yTi-ab — ^^'^x'^ab = £0^^^^)'^^'^ Qab, (4-26) 

where the vector field 1) is defined by 



(7, 



(1.1) = ^^1,1) + ^[^(0,1), ^'^'Ixr + ^[e(i,o)/°''i^]^ (4-27) 



Here again, we have the gauge transformation rule H4-26p . which has the same form 
as that for linear perturbations given in ()4-2() and (|4-3() . Then, as in the previous 

cases of (2, 0) and (0, 2) order perturbations, we can decompose the variables ^^'^Wafe 
as 

^^'^lUab =: ^^'^inab + Vj^'^iXb + Vb^^'^iXa, (4-28) 
where ^^'^^Tiab is gauge invariant and ^^'^^Xa transforms as 

^'''^X- = '•'''ix- + e^'i,!) + l[^^o,i),^'ixr + ^[^(1,0), (4-29) 

4.1.3. Third order metric perturbation 

Finally, we give the definitions of the gauge invariant variables for third order 
metric perturbations. The procedure to define gauge invariant variables is completely 
parallel to that in the case of second order perturbations. First, we define the 
variables ^^''^^'Hab at each order in such a way that their gauge transformation rules 
have the same form as those for the linear perturbations. If we define the variable 
we can extract the gauge invariant part of the higher order perturbation 
theory from this variable using the procedure employed for the linear perturbations 
of the metric. 

The non-trivial point of this procedure lies only in defining the variables ^^'''^Tiab- 
The remaining part is trivial, due to the assumption that there exists a procedure 
for the linear perturbation theory. The definitions of the variables ^^''^^'Hab of O(A^), 
0(A2e), 0(Ae2) and 0{e^) are as follows: 

(3.0) ._(3,0), (2,0), „ /> (1.0), , n n2 (1,0), 
yiT-ab ■— v^ab " •J-*' (1,0) ^ v"-afe " •J-*' (2,0) ^ v"-afe + -J"*- (1,0) yl^ab 

+ |3i^(i,o^^,£(2,o)^^ — °^fi,o)^| ^"■'s'afe, (4-30) 

(2.1) .-(2,1);, 9/> (1,1);, n (2,0), 

x^ x^ 

+ I 2£(l,0) y£(1,1) V + "^(0,1) y£{2,0) y 

V x^ x^ x^ x^ 
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-£(1,0) y£{0,l) y£{l,0) y 



(0) 



(0,3)^ 
X ' '-ab 



x^ 

^ ' X^ab — 2£(o,i) ' xhab — £(1,0) ' x^ab 



dab, 
(0,2) , 



- £ 



-(0,2) 



X 



£ 



(0,1) 



X 



(1,0), 

x' 



2i £ 



(1,1) 

X 



(1,0) v£(0,l) Y 
X^ x^ 



+ £ 



(0,1) 



X 



£ 



(1,0) 



X 



+ i2£ (0,1) y£(l,l) y + £(1,0) ^£(0,2) „ 
x^ x^ x^ x^ 



-£(0,1) £(1,0) £(0,1) ^ 

x^ x^ x^ 



(0) 



(0,3) , 
X'''ab 



q p (0,2) , 

'3 J- (0,1) X^a 



dab, 

(0,: 
X X' 



(4-31) 



(0,1), 
X"'ab 



(4-32) 



3£(o,2) Y^ ' xhab + 3£^o,i) 



X 



(0,1), 
X'^ab 



+ <j 3£(o,i) ^£(0,2) ^ 



£ 



(0,1) 



X 



(0) 



9ab 



(4-33) 



with the gauge X. The variables ^^'"yHab for p + g = 3 are defined for the gauge y 
in same forms as Eqs. H4-3Up - ()4-33p with the replacement X ^ y. Using the gauge 
transformation rules (|T^ - l|nn|l . (jEHl), (|^^^ and (jT^ . we can easily check that 
the variable ^^'^T^afe (p + (7 = 3) is transformed as 



yli-ab X ''O.b — ■^cr, 



(p,9) yabi 



(4-34) 



where the vector fields 0""^ are defined by 



fa 

•^(3,0) ■- ?(3, 



0) 



(T, 



— C« 
(2,1) — ?(2,1) 



+ 



<7, 



(1,2) 



'(1,2) 



+ 



■^(2,0) + ^^''a'-'^' '^(1,0) 



+ 

"(0,1) 



2^(1,0) + 



(1,0) 



+ 



+ 



A'^' ^(1,0) 
(1,0) 



X 



(1,0) Y 
X^^ 



X, C(i,o) 



^(1,0)/"";^-^ 

(2,0) ^ ^ 
- x^ ~ M2,0)>?{0,1) 



^(1,0)' 



(0,1) 
X 



X 



(0,2) Y i: i: 
X^ ~ M0,2):?{1,0) 



+ 
+ 



(0,1) Y C 



(0,1) Y 
X^^ 



?(0,1)> 



(1,0) 



X 



X 



(0,3) •= 



^(0,3) 



' (0,2) Y c 

4(0,2)+ A'^' 4(0,1) 



+ 



24(0,1) 



(0,1) Y 
X^^ 



C(0,1) 



(0,1) 
' X 



X 



(4-35) 



(4-36) 



(4-37) 



(4-38) 



The gauge transformations given in ()4-34() imply that we can decompose the 
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variables ^''^^Hab as 



M-Hab := ^i-Hai, + Va^iXb + V,^'$X„, (4-39) 

using the procedure to find the gauge invariant variables for linear order pertur- 
bations, where ^'^Ha6 is gauge invariant and ^''^X^ at each order is transformed 
as 

a 



(3,0) 



(2,1) 



(3,0) 



^ -^(3,0) 
+ 



< , (2,0) Y c 

?(2,0) + A'-^'Ml.O) 



9^ , (1,0) ^ , (1,0) y 

^qi,0) + A"^' ?(1,0)) 



(4-40) 



(2,1) 



- S(2,l) 
+ 

+ 



+ 



^(1,0) - 2^^'!:^.^ - 2^(1,1) , ^(1,0) 



(1,0) Y f 



^^'5^-^(1,1) '^(0,1) 



(1,0) Y 



qi,o)i A"^ 



(4-41) 



(1,2) 

y-^ 



(1,2) 
AT 



- ^(1,2) + ^(0,1) - 2''''';^X - 2^(1^1) , ^(0,1) 



(1,1) 



+ 

+ 



(0,1) 



^^,C(1,0) - - C(0,2),^(l,0) 



(0,1) Y 



t (1-0) 
?(0,1)> A" 



(0,3) y-a 

y-^ 



(0,3) 



AT^ - ?(0,3) 



3 I6n,2) + ^°'^^,^(o,i) 



+2 



[2^(0, 



1) 



, (0,1) Y 



?(0,1)! A-^ 



(4-42) 



(4-43) 



Thus, we have found recursively the gauge invariant variables for higher order 
metric perturbations up to third order. 

4.2. Gauge invariant variables for matter perturbations 

As shown above, we can find the gauge invariant variables of higher order metric 
perturbations. To do so, we have decomposed the metric perturbations at each order 
into the gauge invariant variables ^^'^Haft and the gauge variant vector variables 

^^''^■^Xa- The gauge variant parts, ''^'^X^ of the 0{X^e'^) metric perturbations are 
irrelevant as physical metric perturbations. However, using these gauge variant 
parts, we can define the gauge invariant variables for the physical fields, other than 
the metric. 

Here, we give explicit definitions of gauge invariant variables for perturbations 
of an arbitrary tensor field Q up to third order: 



(4-44) 
(4-45) 



4''°^ 2 ■■= ^x'^Q - 2^(1.0) ^4'°^Q - \ ^i-:") X - "^(1.0) ^ \ Qo, (4-46) 
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^^^^\x ~ 2'^*^'"^x°^<"'^j.x ~ '2^^^-\x^^^-°\x^ (4-47) 

40'2)q = _ 2£,„_,, ^^O'Dq _ 1^(0.2) ^ - "^^?o,2) ^} go, (4-48) 

(^^'°^Q := b^^^^Q — 3£(i,o) y^b^^f'^Q — 3 |<£(2,o) ^ — i^fi.o) ^| b^x^^ Q 

+ 4,,. Iqo, (4-49) 



x^ 



-2 .£(1,1) ^ — -£{ifi) y£{o,i)^^ — -£(o,i^ji^<£(i,o^^ I b^^''^^Q 



\ £(2,1) y — 2£(l,0) £(1,1) — <£(0,1) ^=£(2,0) V 

+£(1,0) ^£(0,1) ^£(1,0) ^ }Qo, (4-50) 



<5v'^^Q :— b^^''^^Q — 2£(o,i) (5^'^''Q — ■^(i.o) yb^v''^^Q 

x^ x^ 

— { £(0.2) Y — £%,!) ^1 b^^'^^'Q 

1 „ . 1 



-2 <i £(1,1) ^ - -£(i,o^^£(o,i^^ - -£(o,i^^£(i,o^^ j> 4°'^^g 



.£(1,2) „ — 2£(0.1) ^£(1,1) y — £(1,0) y£(0,2) „ 
x-^ X^ X^ X^ X^ 



-{■ 

+ £(0,1) ^£(1,0) ^£(0,1) ^1 Qoi (4'51) 
b^x^^'Q := b^^'^^Q — 3£(o,i) — 3 < £(o,2) ^ — £%,i) ^ \ b^^'^^Q 



£(0,3) y — 3£(o,i) y£(o,2) y + £(0,1) ^ > Qo- (4-52) 

x^ x^ x^ x^ 



Straightforward calculations using the gauge transformation rules p-30|) - H3-38p . (|4-13|) . 



(ITOl) . (j3211), (irai) and show that these variables jJ'^^Q are gauge 



invariant. 

§5. Summary and Discussions 

In this paper, we have presented the procedure to find gauge invariant variables 
of two-parameter nonlinear metric and matter perturbations. To describe this proce- 
dure, we have assumed the existence of the analogous procedure for linear perturba- 
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tions. As emphasized in the main text, the decomposition of the Unear perturbation 
of the metric into gauge invariant and gauge variant parts is non-trivial. It depends 
crucially on the background spacetime. However, if we assume the existence of this 
procedure, we can always define the gauge invariant variables for higher order metric 
and matter perturbations. 

The procedure presented above is summarized as follows. Suppose that gauge 
transformation rules for a (p, g)-order metric perturbation ^^''^^hab are given by 

^^''yhab — ^^''^^hab = ^[^{p,q)^ ^ '''l;hab], (5-1) 

where F is a function determined by the gauge transformation rule Af — > 3^ of the 
(p, g)-order metric perturbation ^^''^^hab, and i and j are integers that satisfy the 
conditions i < p, j < q and i + j ^ p -\- q. Further, suppose that for any i and j, we 
have already defined the tensor fields ^^''^j^Ti.ab and the vector fields ^^'"'^X" as in the 
main text. Then, defining the variable ^^''^j^Tiab for the 0{\^e'^) metric perturbation 

by 

there exists a vector field 0""^ such that the variable ^^'''■^'Hab transforms as 

y'Hab ^ x'Hab — "^""(p.q)^ ^ 9ab (^'3) 

under the gauge transformation X ^ y. Then, using the same procedure as for 
linear perturbations, we can decompose ^^'"^^Hab as 

Xl^ab—- x'^ab + ^a x^b+^/b X^a, 

where the variables ^^'"^^Tiab and ^^''^Xq are the gauge invariant and gauge variant 
parts of the 0(A^e'') metric perturbation, respectively. Because the gauge transfor- 
mation rule for the matter perturbations 5^''^^Q is given by 

4^'^)g - 6^r^Q = (5-5) 

the corresponding gauge invariant variables are defined by 

55^''?)q = 5j''?)Q + -^'''lx%6^i''^Q]. (5-6) 

A procedure similar to that presented here for second order perturbations in 
the one-parameter case was previously obtained by Campanelli and Lousto.l^ They 
applied that procedure to second order perturbations of a Kerr black hole. We have 
confirmed their procedure up to third order in the two-parameter case. Though the 
gauge transformation rule for the perturbations of arbitrary order is not yet known, 
we conjecture that the procedure considered here is applicable to arbitrary order 
perturbations. We also believe that this procedure can be confirmed by induction, 
once the gauge transformation rule for an arbitrary order is obtained. We leave this 
for a future work. 
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In addition to the interesting mathematical framework, it is also interesting to 
apply it to the astrophysical systems, such as oscillating relativistic rotating stars. 
Many astrophysical systems can be described well by perturbation theory with two 
parameters. One merit of applying the procedure presented here to such systems 
is the gauge ambiguities are removed. When we apply this procedure to oscillating 
stars, a careful analysis is necessary to properly treat the boundary conditions at 
the surface of the star and the displacement of this surface. Similar situations are 
considered in preveous papers coauthored by the present author.*^ In those papers, 
a particular gauge fixing is necessary to match the perturbative solutions at the 
boundary of the surface of the matter distribution when we construct global solutions 
and when we define the perturbative displacement of the matter surface. Similar 
problems also arise in the perturbative analysis of spacetimes with boundaries, such 
as brane worlds.!^ In the investigation of these situations, the gauge transformation 
rules derived in BGS2003 and the gauge invariant variables defined here become 
powerful tools. Because of their applicability to various situations, we also expect 
that the techniques developed here will play a key role in progress of theoretical 
physics. 
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Appendix A 

Equivalence with the Representation of Bruni et al. 

In this appendix, we derive the representation of the coefficients of the formal 
Taylor expansion ()2-3() of the pull-back of a diffeomorphism in terms of the suitable 
derivative operators ^(p^q)- As mentioned in the main text, for each of these opera- 
tors, there is a vector field such that the operators satisfy Eq. (|2-13j) . The existence 
of these vector fields is guaranteed by the fact that the operator >C(p^g) of each order 
satisfies the Leibniz rule. From this fact, the definitions of the derivative operators 
H2-4() - H2-12() are obtained. As shown in the Appendix in BGS2003, the representa- 
tion of the Taylor expansion of ^\^f for an arbitrary function / can be extended 
to that for an arbitrary tensor field Q on Ad. Therefore, in this appendix, we only 
consider the Taylor expansion of the pull-back ^\^f for an arbitrary scalar function 



fe,fc'=0 I. ) A=e=0 



where J-{M) denotes the algebra of functions on M.. 

Although the operators d/dX and d/de in the bracket {*}A=e=o of Eq. HA-1|) are 
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simply symbolic notation, we stipulate the properties 
gn+l ^ f 9 / 9" 



^" ^ ^-^le/ ) ^ , (A-2) 

A=e=0 



dX^ \dX 



8,„,8/Wn . (A3) 



A=e=0 



A=e=0 

^ ^ ^ '^I.e/) I , (A-4) 



de \dX 



for V/ G J-{Ai)^ where n is an arbitrary finite integer. These properties imply that 
the operators d/dX and d/de are in fact not simply symbolic notation but indeed 
the usual partial differential operators on R,^. We note that the properties (|A-5|) 
and (|A-6|) are the Leibniz rules, which play important roles when we derive the 
representation of the Taylor expansion (jA-ip in terms of suitable Lie derivatives. 

We can easily show that the derivative operators >C(i o) ^-^id '^(0,1)1 defined by 
(|2-4j) and (|2-5|) . respectively, satisfy the Leibniz rule 

= 2//:(p,,)/ (A-7) 



due to the properties ()A-5|) and ()A-6|) . In the higher order coefficients of the expan- 
sion (|AT|) . the properties (|A-5|) . (|A-6|) and (|A-7|) lead to non-trivial combinations of 
the linear operators. In BGS2003, it is commented that the representations of the 
higher order coefficients is not unique, and the following combinations are derived: 

§X^lj] =^imf, (A-8) 



A=e=0 



-^A f\ 



A=e=0 



^(0,1)/, (A-9) 

{^^A,e/}^^^^^ = ^(1,1)/ + (eo/:(l,0)^(0,l) + 61^(0,1)^(1,0)) /, (A-11) 



A=e=0 

'^*xJ\ = Ao,2)/ + 4o,i)/' (A-12) 



92 

-j^'^X,eJ ( - H0,2)J ^^(0,1) 

^6 ' A=e=0 
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^(3,0)/ + 3/^(i,o)'C(2,o)/ + ^(1,0)/' 

^(2,1)/ + 2i2(l,o)'C(l,l)/ 

+'^{o,i)'^(2,o)/ + 2e2'C(i,o)'C(o,i)>C(i,o)/ 
+(ei - e2)/:(o,i)4i,o)/ + (eo - €2)4,0)^(0,1)/, (A-14) 

^(l,2)f + 2/^(0,l)'C(i,i)/ 

+'C(i,o)'^(o,2)/ + 2e3'C(o,i)'C(i,o)'C(o,i)/ 
+ (eo - e3)/:(i,o)4,i)/ + («i - «3)'C(o,i)'^(i-0)/> (A-15) 

^(0,3)/ + 3£(o,i)/:(o,2)/ + 4o,i)/- (A-16) 

Here, the parameters (i = 0, 1, 2, 3) are constants that satisfy the condition 

eo + ei = l. (A-17) 

These parameters result from the fact that the representation of the higher order 
coefficients in terms of the derivative operators is not unique. 

As emphasized in BGS2003, the parameters ej in the representations (|A-8|) - 
(|AT6|) are not essential. Each operator >C(p,g) is regarded as the Lie derivative with 
respect to the vector field ^-^ as mentioned above (see Eq. (|2-13() '). Once we obtain 
the representation of ^[p^q) in terms of the Lie derivative, the parameters are 
always eliminated through the replacement of the vector field. As a result, we obtain 
the ^''canonical representation" ()2-14|) in the main text. Therefore, the "canonical 
representation" H2-14() of the Taylor expansion with two infinitesimal parameters is 
equivalent to that in BGS2003 as shown below. 

To show this, it is necessary to give an explicit derivation of the representations 
(|A-8|) ~ HA-16|) . The derivations of these representations are done recursively from 
lower order representations. Because the derivation of each representation is similar, 
it is enough to present the explicit derivation of the representation of the coefficients 
of O(A^e), and we start from the point where the lower order representations HA-8() - 
(|AT3|) are already given. It is obvious that the coefficient of O(A^e) in the expansion 
HAT|) includes a linear combination of the following terms: 

'C(0,l)^(2,0)/> ^{2,0)^(0,1) f 7 ^(1,0)^(1,1) fi ^{1,1)^(1,0) 

-^(0, 1)^(1,0)^(1,0)/) 'C(i,o)'C(o,i)'C(i,o)/) ^(1,0)^(1,0)^(0,1) f- (A-18) 
Then, we start from the following form of the derivative operator 'C(2,i)- 

^(^'^)^^={5^^^'^4a=.=o 

+Pl'C(0,l)'C(2,0)/ +^2-^(1, 0)^(1,1)/ 

+P3^(0,l)^(l,0)^(l,0)f + 1'4^(1,0)^(0,1)'^(1,0)/ + ^5-^(1, 0)^(1, 0)^(0,1)/ 

+Pe^(2,o)^(o,i)f + P7^(i,i)^(i,o)f- (A-19) 



l^n f] = 
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Following the argument in BGS2003, we choose 

P<i=P7 = 0. (A-20) 

This choice is always possible. The lower order derivative operators ^(p^q) are re- 
garded as the Lie derivative with appropriate generators as in Eq. (|2 13j) and 
the first term on the second line and the first term on the fourth line in the right-hand 
side of Eq. (|A-19|) are given by 

= (pi+P6)i?C(o,i)-^^C(2,o)/ + -fp6K(2,o),5(o,i)]/- (^-21) 

Because, as we will show, jC(2,i) = £^1^2 1) ^ven if 7^ 0; we can always choose Pq = 
by making the replacement 

Pl+P6^Pl, C('2,l)+P6[?(2,0),C{0,l)]"^e('2,l). (A-22) 

without loss of generality. The same argument can be applied to the parameter p^. 
Therefore, we may choose Pq = Pi = 0, without loss of generality. We note that 
similar arguments can be applied to the cases of the parameter ^3, p4 and p^. We 
can fix these parameters by using the replacement of the vector field i)- However, 
we do not proceed with this argument here, because we have confirmed that the 
representation in BGS2003, i.e. (|X^ - (|XT6)) in this paper, and Eqs. ^T^-^TV^ 
are equivalent, at least up to fourth order. 

To guarantee the Leibniz rule for the derivative operator C(2^i), we first consider 
the operation 

Using Eqs. we obtain 

,2 



A=e=0 



e 



Substituting Eqs. ^iKM - ^^KA% . (|X49l) and into Eq. (IX^ . we obtain the 

representation of Eq. HA-23|) . which includes the term 2/£(2,i)/- Oii the other hand, 
a direct calculation from Eq. (|AT9|) with Eqs. ()A-2fl|) gives the representation of 
Eq. HA-23|) which includes the term >C(2,i)/^- To obtain this second representation 
of Eq. ()A-23|) . we can use the Leibniz rules for the derivative operators, except for 
£(2,1)) because the Leibniz rules for the lower order operators can be demonstrated 
by arguments similar to that given here, and these can be confirmed before applying 
the arguments for the operator C(2,i)- From these two representations of Eq. ()A-23|) . 
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we obtain 

J0i2,i){f) - 2/£(2,i)/ = 2(m + i)(/:(o,i)/)(^(2,o)/) 

+2(p2 + 2)(/:(i,o)/)(/:(i,i)/) 

+ 2(P3 +P4 +P5 + i)(/:(o,i)/)(4,o)/) 

+2(p4 + 2p3 + 2ei)(/:(o,i)/:(i,o)/)(>C(i,o)/) 

+2(p4 + 2p5 + 2eo)(/:(i,o)/:(o,i)/)(>C(i,o)/). (A-25) 

This shows that the derivative operator >C(2,i) satisfies the Leibniz rule jC(2,i)(/^) = 
2/^(2,1)/ an arbitrary function / G ^(-M) iff 

Pi = -1, P2 = -2, P3 = -(ei-e2), P5 = -(eo-e2), (A-26) 

where eo and ei appear from the representation of through the substitution of 

Eq. ()A-11|) into Eq. (|A-24|) . and €2 is an arbitrary parameter that is not determined 
by the Leibniz rule for ^(2,1)- Thus, the representation HA-14|) is obtained. 

Finally, we show that the undetermined parameter {i = 0, • • • , 3) can be 
eliminated through the replacement of the vector field s ^ -^ . First, we consider the 
derivative operator HA-11|1 of 0(Ae). Then, using Eq. (|AT7|) . the replacement 

^(1,1) = + ^(eo - el)[e(l,o),e(o,l)]^ (A-27) 

leads 

{ 5Si^^'^4A=.=o ^ ^''^^^^^^^ (^^«(^.o)^«(o.i) + ^^«(o.) A(i,o)) /• (A-28) 
This implies that the parameters eo and ei should be chosen so that 

eo = ei = ^ (A-29) 
without loss of generality. Similarly, the replacements of the generators ^(2,1) and 

C(l,2), 

<2,i) ■■= ^(2,1) + (ei - e2)[e(i,o), [e(l,o),e(o,l)]]^ (A-30) 

^(1,2) •= ^(1,2) + («i - ^3) [^(0,1) ' K(o,i) > ^{1,0)]]": (A-31) 

lead to 

"'"°^C(o,i) °^?(2,o) + °^5(i,o)'^C(o,i)'^5{i,o)' (A-32) 

"'"°^?(i,o) °^C(o,2) + ■^5(0,1) "^^ci.oj-^Sco,!)' (A-33) 



Gauge Invariant Variables in Two- Parameter Nonlinear Perturbations 33 



These show that the parameters 62 and €3 should be chosen so that 

£2 = £3 = ^ (A-34) 

without loss of generality. Thus, the undetermined parameters in the representation 
derived in BGS2003 are eliminated through the replacement of generators ^(^q)) and 
the representation of the Taylor expansion in BGS2003 is equivalent to (|2-14|) in the 
main text. This also implies that the representation of the Taylor expansion of the 
pull-back is not unique, but this non-uniqueness causes no serious problems. 

References 

1) D. Kramer, H. Stephani, M. A. H. MacCallum and E. Herlt, Exact Solutions of Einstein's 
Field Equations (Cambridge: Cambridge University Press, 1980). 

2) J. M. Bardeen, Phys. Rev. D 22 (1980), 1882. 

3) H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. No.78 (1984), 1. 

4) V. F. Mukhanov, H. A. Feildman, and R. H. Brandenberger, Phys. Rep. 215 (1992), 203. 

5) S. Chandrasekhar, The mathematical theory of black holes (Oxford: Clarendon Press, 
1983). 

6) R. J. Gleiser, C. O. Nicasio, R. H. Price and J. Pullin, Phys. Rep. 325 (2000), 41. 

7) K. D. Kokkotas and B. G. Schmidt, Living Rev. Relativity 2 (1999), 2. 

8) Y. Kojima, Prog. Theor. Phys. Suppl. No.l28 (1997), 251. 

9) J. B. Hartile, Astrophys. J. 150 (1967), 1005. 

S. Chandrasekhar and J. C. Miller, Mon. Not. R. Astron. Soc. 167 (1974), 63. 

10) N. Stergioulas, Living Rev. Relativity" 2, (1999) 2. 

11) R. M. Wald, General Relativity (Chicago, IL: University of Chicago Press, 1984). 

12) J. M. Stewart and M. Walker, Proc. R. Soc. London A 341 (1974), 49. 
J. M. Stewart, Class. Quantum Grav. 7 (1990), 1169. 

J. M. Stewart, Advanced General Relativity (Cambridge University Press, Cambridge, 
1991). 

13) K. Nakamura, Prog. Theor. Phys. 110, (2003), 201. 

14) M. Bruni, L. Gualtieri, and C. F. Sopuerta, Class. Quantum Grav. 20, (2003), 535. 

15) M. Bruni, S. Matarrese, S. MoUerach and S. Sqonego, Class. Quantum Grav. 14, (1997), 
2585. 

16) Y. Choquet-Bruhat, C. DeWitt-Morette and M. Dillard-Bleick, Analysis, Manifolds and 
Physics (North-Holland, Amsterdam, 1982). 

17) S. Sonego and M. Bruni, Commun. Math. Phys. 193, (1998), 209. 

18) W. Thirring, A Course in Mathematical Physics: I. Classical Dynamical Systems (New 
York: Springer, 1978). 

19) S. Kobayashi and K. Nomizu, Foundation of Differential Geometry vol. I (Wiley, New 
York, 1963). 

20) S. Matarrese, S. MoUerach and M. Bruni, Phys. Rev. D58, (1998), 043504. 

21) M. Bruni and S. Sonego, Class. Quantum Grav. 16, (1999), L29. 

22) U. H. Gerlach and U. K. Sengupta, Phys. Rev. D 19, (1979), 2268. 
U. H. Gerlach and U. K. Sengupta, Phys. Rev. D 20, (1979), 3009. 
U. H. Gerlach and U. K. Sengupta, Phys. Rev. D 22, (1980), 1300. 
U. H. Gerlach and U. K. Sengupta, J. Math. Phys. 20, (1979), 2540. 

23) M. CampaneUi and C. O. Lousto, Phys. Rev. D 59, (1999), 124022. 

24) K. Nakamura, A. Ishibashi and H. Ishihara, Phys. Rev. D62 (2000), 101502(R). 
K. Nakamura and H. Ishihara, Phys. Rev. D 63 (2001), 127501. 

25) L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999), 4690. 
L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999), 3370. 
J. Garriga and T. Tanaka, Phys. Rev. Lett. 84 (2000), 2778. 



